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ABSTRACT 

We recalculate the diffuse ffux of high energy neutrinos produced by Gamma 
Ray Bursts (GRB) in the relativistic fireball model. Although we confirm that 
the average single burst produces only ~ 10~ 2 high energy neutrino events in 
a detector with 1 km 2 effective area, i.e. about 10 events per year, we show 
that the observed rate is dominated by burst-to-burst fluctuations which are 
very large. We find event rates that are expected to be larger by one order of 
magnitude, likely more, which are dominated by a few very bright bursts. This 
greatly simplifies their detection. 

Subject headings: gamma rays: bursts 
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1. High Energy Neutrinos from Relativistic Fireballs 

The evidence has been steadily accumulating that GRB emission is the result of a 
relativistically expanding fireball energized by a process involving neutron stars or black 
holes ( Piran 1999|) . In the early stage, the fireball cannot emit efficiently because the 
radiation is trapped due to the very large optical depth. The fireball's energy is dissipated 
in kinetic energy until it becomes optically thin and produces the observed display. This 
scenario can accommodate the observed energy- and time scales provided the bulk Lorentz 
factor 7 of the expanding shock is ~ 300. 

The production of high energy neutrinos is anticipated: protons, accelerated in the 
kinetic phase of the shock, interact with photons producing charged pions which are the 
parents of neutrinos ( |Waxman fc Bahcall 199TD - Standard particle physics and fireball 
phenomenology are sufficient to compute the neutrino flux ( |Waxman fc Bahcall 1997| ; 
[Halzen 1998| ) as well as the observed rates in high energy neutrino telescopes. 

The observation of GRB neutrinos over a cosmological baseline has scientific potential 
beyond testing the "best-buy" fireball model: the observations can test with unmatched 
precision special relativity and the equivalence principle, and study oscillating neutrino 
flavors over the ultimate baseline of z ^ 1. 

The anticipated neutrino flux traces the broken power-law spectrum observed for the 
photons, which provide the target material for neutrino production: 

= -A- ioiE<E B (1) 

4> = A for E>E B , (2) 

where A is a normalization constant which is determined from energy considerations and 
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Eb ~ 700 TeV ( Halzen 1998| ). The total energy in GRB neutrinos is given by 



c 1 

^tot = 1 — Tc/tt bubble PE i (3) 
47T Z 

where f n is the fraction of proton energy that goes into pion production, ihubbie is 10 Byrs 
and pE is the injection rate into protons accelerated to high energies in the kinetic fireball. 
This is a critical parameter in the calculations; it can be fixed, for instance, by assuming 
that GRB are the source of the highest energy cosmic rays beyond the ankle flWaxman 



|1995; Milgrom fc Usov 1995; Vietri 1995|) in the spectrum, or p E ~ 4 x 10 44 erg Mpc" 3 yr _1 
( [Halzen 1998|) . The factor f n represent the fraction of total energy going into pion 
production. It is calculated by known particle physics and is of order 15% ([Waxman fc 



Pahcall 1997] ; [Halzen 1998|) . We remind the reader that these assumptions reproduce 



the observed average photon energy per burst if equal energy goes into the hadronic and 
electromagnetic component of the fireball ( |Waxman fc Bahcall 1997 ). 



Now, normalizing the flux for this total energy: 

A f E B fEm^ dE , . 

Fm = ~ET dE + A — . (4) 

&B JE min JE B Hi 

Approximating E B — E min ~ E B , the integration constant is found to be 
1.20 x 10~ 12 TeVcm -2 s _1 sr _1 . The quantities /„- and E B are calculated using Eqs. (4) 
and (5) in Waxman and Bahcall (1997). The GRB parameters entering these equations 
are chosen following Halzen (1998) and assumed to be independent of 7. We assumed 
£ max = 10 7 TeV. 

Given the neutrino flux (J), the event rates in a neutrino detector are obtained by a 



straightforward method QHalzen 1998 



N 



J tf>{E)P v ^dE , (5) 



where Pp^^ is the detection probability for neutrinos of muon flavor. It is determined, as 
a function of energy, by the neutrino cross sections and the range of the secondary muon 
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( paisser, Halzen fc Stanev 1995; Gandhi et al. 1996 ). We obtain an observed neutrino flux 



of order 10 yr 1 km 2 assuming 10 3 bursts per year. This result is somewhat lower but not 
inconsistent with those obtained previously. 

One should keep in mind that neutrinos can also be produced, possibly with higher 
fluxes, in other stages of the fireball, e.g. when it expands through the interstellar medium 
( |Katz 1994; Halzen fc Jaczko 1996; Vietri 1998; Boettcher fc Dermer 199SQ . 



2. Burst-to-Burst Fluctuations 

We here want to draw attention to the fact that this result should not be used without 
consideration of the burst-to-burst fluctuations which are very large indeed; see also 
Permer, Boettcher fc Chiang (1999)| . We will, in fact, conclude that from an observational 



point of view the relevant rates are determined by the fluctuations, and not by the average 
burst. Our calculations are performed by fluctuating individual bursts according to the 
model described above with i) the square of the distance which is assumed to follow a 
Euclidean or cosmological distribution, ii) with the energy assuming that the neutrino rate 
depends linearly on energy and follows a simple step function with ten percent of GRB 
producing more energy than average by a factor of ten, and one percent by a factor of 100, 
and, most importantly, iii) fluctuations in the 7 factor around its average value of 300. 
The fluctuations in 7 affect the value of /„- which varies approximately as 7~ 4 ([Waxman & 



Bahcall 1997| ) as well as the position of the break energy which varies as 7 2 . For a detailed 



discussion see Halzen (1998). Both effects are taken into account. Clearly a factor of ten 
variation of 7 leads to a change in flux by roughly 4 orders of magnitude. 

The origin of the large fluctuations with 7 should be a general property of boosted 
accelerators. With high luminosities emitted over short times, the large photon density 
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renders the GRB opaque unless 7 is very large. Only transparent sources with large boost 
factors emit photons. They are however relatively weak neutrino sources because the actual 
photon target density in the fireball is diluted by the large Lorentz factor. 

This raises the unanswered question whether there are bursts with lower 7-factors. 
Because of absorption such sources would emit less photons of lower energy and could 
have been missed by selection effects; they would be spectacular neutrino sources. Some 
have argued ( [Stern 19991 ) f° r (7) ~ 30 on the basis that the unusual fluctuations in the 
morphology of GRB can only be produced in a relatively dense, turbulent medium. 

3. Monte Carlo Simulation Results 

We will illustrate the effect of fluctuations by Monte Carlo simulation of GRB. For 
a Euclidean distribution the calculation can be performed analytically; in other cases 
the Monte Carlo method evaluates the integrals. The overwhelming effect of fluctuations 
is demonstrated in sample results shown in Figs. 1-3 and Table l,a-c. Not knowing the 
distribution of 7-factors around its average, we have parametrized it as a Gaussian with 
widths a, a' below and above the average value of 300. We chose a' to be either or 300, 
to illustrate the effect of allowing GRB with Lorentz factors up to 10 3 with a Gaussian 
probability. The critical, and unknown, parameter is a. It may, in fact, be more important 
than any other parameter entering the calculation. As far as we know, neither theory nor 
experiment provide compelling information at this point. Note that we require a > 70 in 
order to allow a significant part of the bursts to have 7-factors less than 10 2 , or one third the 
average value; see Table 1. The value of a' is less critical. A value of 300 allows for Lorentz 
factors as large as 10 3 . The dominant effect of this is a renormalization of the neutrino 
rates because a fraction of the bursts now have a large 7-factor and, as a consequence, a 
low neutrino flux. We discuss some quantitative examples next. 
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Firstly, even in the absence of the dominant fluctuations in 7, the rate of 9 detected 
neutrinos per year over 10 3 bursts, becomes 30-90 in the presence of fluctuations in energy 
and distance. The range covers a variety of assumptions for GRB distributions, which 
range from Euclidean to cosmological; the issue is not important here because other factors 
dominate the fluctuations. Every 2 years there will be an event with 7 neutrinos in a single 
burst. For a = 70, the rate is ~ 600 per year in a kilometer square detector, with 23 
individual bursts yielding more than 4 neutrinos and 7 yielding more than 10 in a single 
year! The results for other values of a are tabulated in Table 1. The number of events per 
year for a range of values of a is shown in Fig. 1. The frequency of GRB producing seven 
or more neutrinos is shown in Fig. 2. The neutrino multiplicity of bursts for a = 30, 60 
and 75 is shown in Fig. 3. Note that absorption in the source eventually reduces the overall 
rates when 7 is much smaller than average, or a larger than 60. This factor is included in 
our calculations ([Halzen 1998| ). 

To this point we have assumed that all quantities fluctuate independently. Interestingly, 
the value of a = 70 is obtained in the scenario where instead fluctuations in 7 are a 
consequence of fluctuations in energy. In order not to double count, fluctuations in energy 
itself, which only contribute a factor of 3 to the neutrino rate anyway, should be omitted. 

Even though bursts with somewhat lower 7 produce neutrinos with reduced energy, 
compared to 700 TeV, over a longer time-scale than 1 second, the signatures of these models 
are spectacular. Existing detectors with effective area at the 1-10% of a kilometer squared, 
should produce results relevant to the open question on the distribution of bulk Lorentz 
factors in the fireball model. 

Independent of the numerics, the fact that a single GRB with high energy, close 
proximity and a relatively low Lorentz factor can reasonably produce more detectable 
neutrino events than all other GRB over several years time, renders the result of the 
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straight forward diffuse flux calculation observationally misleading. Our calculations 
suggest that it is far more likely that neutrino detectors detect one GRB with favorable 
characteristics, than hundreds with average values. Clearly, our observations are relevant 
for other GRB models, as well as for blazars and any other boosted sources. They are also 
applicable to photons and may represent the underlying mechanism for the fact that the 
TeV extra-galactic sky consists of a few bursting sources. We expect no direct correlation 
between neutrino and high energy gamma sources because cosmic events with abundant 
neutrino production are almost certainly opaque to high energy photons. 

4. Conclusions 

Although the average event rates predicted with typical GRB parameters appear 
somewhat discouraging to present and future Cerenkov neutrino detector experiments, 
the fluctuations in these calculations are more significant and affect the prospects for 
detection. We speculated on the distribution of the parameters entering the fireball model 
calculation and used a Monte Carlo simulation to estimate the actual event rates. The 
result of these simulations show that a kilometer scale detector could be expected to observe 
tens or hundreds of events per year. To improve on the reliability of this estimate, a well 
defined distribution for the Lorentz factor must be determined. Contrariwise, not observing 
neutrino GRB after years of observation will result in a strong limit on the number of 
accelerated protons in the kinetic phase of the burst, or in fine-tuned, high values of the 
Lorentz boost factor. 
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Fig. 1. — Rate of high energy neutrino events in a detector with 1km 2 effective area as a 
function of the range of the bulk Lorentz factor 7. The range is assumed to be Gaussian 
with width a below the average of 300, and o' above. The latter is taken to be either 0, or 
300. 
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Fig. 2. — Yearly rate of individual neutrino bursts with more than 7 events in a detector 
with 1 km 2 effective area as a function of the range of the bulk Lorentz factor 7. The range 
is assumed to be Gaussian with width a below the average of 300, and a 1 above. The latter 
is taken to be either 0, or 300. 
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Fig. 3. — Neutrino multiplicity distribution in individual GRB in a detector with 1km 2 
effective area. The range of the bulk Lorentz factor 7 is assumed to be Gaussian with width 
a = 30, 60 and 75 below the average of 300, and a' taken to be 0. 
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Table 1: High energy neutrino events in a detector with 1 km 2 effective area as a function of 
the range of the bulk Lorentz factor 7. The range is assumed to be Gaussian with width a 
below the average of 300, and a' above. The latter is taken to be either 0, or 300. Shown as 
a function of a is the total number of neutrinos, as well as the yearly number of individual 
bursts with more than 4, 7, and 11 high energy neutrinos. For comparison, the middle line 
shows the result when all fluctuations are removed. 

a' = 



Percent of Lorentz Factors by Range GRB > X Events/yr 

cr 0-10 10-50 50-100 100-200 200-300 Events/yr > 4 > 7 > 11 



80 


0.04 


0.15 


1.09 


19.83 


78.9 
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33 


18 


12 
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10 
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Result if energy and distance are not allowed to vary: 
100 8.6 

cr' = 300 

Percent of Lorentz Factors by Range GRB > X Evcnts/yr 

cr 0-10 10-50 50-100 100-200 200-300 Evcnts/yr > 4 > 7 > 11 
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0.01 


0.05 


0.29 


5.3 
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